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CHAPTER  I. 


MICROSTRIP  ANTENNA  STRUCTURE 

Of  the  various  types  of  electromagnetically  coupled  feed  structures, 
the  coplanar  ones,  commonly  referred  to  in  the  literature  as  slot  line  and 
coplanar  waveguide,  are  of  special  interest  because  of  the  ease  with  which 
active  devices  with  a  planar  architecture  are  integrated  with  the  antenna 
elements.  Engineering  efforts  made  during  the  last  period  were  therefore 
concentrated  on  analysing  and  understanding  the  coplanar  structures. 

Our  preliminary  experiments  showed  that  it  was  possible  to  efficiently 
excite  the  microstrip  antenna  with  the  coplanar  waveguide  feed.  Subsequent 
analysis  and  experiments  were  performed  to  understand  the  excitation  mechanism 
with  the  goal  of  developing  an  analytic  model  that  can  predict  the  performance 
of  the  antenna. 

To  simplify  the  analysis,  several  assumptions  are  made.  The  dimensions 
of  the  antenna  are  such  that  the  cavity  model  for  the  -ntenna  can  be 
invoked.  The  slots  in  the  coplanar  waveguide  are  narrow  compared  to  the 
width  of  the  center  conductor.  This  results  in  a  low  impedance  transmission 
line  with  the  electric  fields  concentrated  within  the  slots.  The  coplanar 
waveguide  can  be  analysed  as  two  parallel  coupled  slot  lines  on  which  an 
odd  mode  has  been  excited.  Huygens'  principle  can  be  applied  to  the  slot 
and  one  obtains  an  equivalent  magnetic  current  source,  equal  in  magnitude 
but  opposite  in  phase,  for  each  slot.  If  an  even  mode  has  also  been  excited 
on  the  coplanar  waveguide  the  current  sources  will  not  be  equal  in  magnitude. 

A  basic  insight  into  the  coupling  mechanism  for  a  coplanar  feed  is 
given  by  the  energy  relation  for  electromagnetic  fields  stated  below. 


2 


-  J(H-K  +  E*J)«dS  =  |^J(±ToE2+i-n0H2)dV+ J  oE2dV  +  ^£xF$*d§ 


P$  =  Power  Supplied  to  Antenna 


J. 

f 

l 


(TEoE2+  T^qH2)  dV  =  Total  Stored  Energy 


<yE  dV  =  Total  Dissipated  Power 


£  x  •  d§  =  Total  Radiated  Power 


There  aust  be  a  source  coaponent  in  parallel  with  the  field  component 
that  is  to  be  excited.  The  larger  the  aean  value  of  the  scalar  product 
of  the  source  and  field  components  is  over  a  closed  surface,  the  greater 
the  power  supplied  to  the  region  under  consideration.  For  example,  by 
feeding  the  patch  antenna  by  a  coplanar  waveguide  as  indicated  in  Fig.  1. 
it  is  possible  to  excite  the  (O.i)  mode.  In  the  (0,1)  mode  the  magnetic 
field,  H.  is  oriented  in  the  x  direction,  and  it  goes  to  zero  on  the 
magnetic  wails  at  y*0  and  y«b.  It  has  maximum  magnitude  at  y=b/2.  There 
is  no  variation  in  the  x  direction  inside  the  cavity  for  the  (0,1)  mode. 
The  electric  fields  are  in  the  z  direction  and  are  zero  along  the  line 
y*b/2  inside  the  cavity.  At  the  end  of  the  coplanar  waveguide  feed,  the 
equivalent  magnetic  current  source,  K,  is  oriented  in  the  x  direction. 

This  source  component  is  in  parallel  with  the  magnetic  field.  H,  of  the 
(0.1)  mode  and  should,  therefore,  excite  that  cavity  mode.  The  y  directed 
components  of  the  source.  K.  should  have  little  effect  on  the  dominant 
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(0,1)  aode  but  may  excite  other  modes.  Maximum  excitation  should  occur 
when  the  feed  is  terminated  at  the  center  of  the  patch  where  the  magnetic 

field  is  maximum.  As  the  feed  is  terminated  closer  to  the  edges  of  the 

patch,  the  excitation  should  go  to  zero  as  the  magnetic  field  goes  to 
zero.  Note  that  excitation  of  the  cavity  occurs  along  the  line  y=b/2 
even  though  the  electric  field  is  zero  along  that  line.  Table  1  lists 
antenna  input  impedance  versus  feed  location.  These  impedance  values 
have  been  measured  at  a  plane  extended  several  wavelengths  from  the  edge 
of  the  patch  antenna.  At  this  point,  to  build  a  good  qualitative  picture 

of  input  impedance  versus  feed  location,  the  magnitude  of  the  input 

impedance  is  the  parameter  of  primary  interest  and  the  phase  shift  due  to 
the  extended  measurement  plane  is  of  secondary  importance.  (Time  domain 
ref lectometry  was  used  to  insure  a  near  perfect  transition  from  a  coaxial 
transmission  line  to  the  coplanar  waveguide  feed,  so  that  the  measured 
impedance  differs  from  the  actual  antenna  impedance  by  only  a  phase  shift 
factor .  ) 

As  one  would  expect,  the  data  show  that  much  of  the  incident  power 
is  reflected  as  the  feed  termination  is  moved  to  the  edge  of  the  patch 
in  the  y  direction.  There  is  little  variation  as  the  feed  is  moved  along 
the  x  direction.  Maximum  power  transfer  from  a  50  ohm  system  to  the 
antenna  occurs  when  the  feed  is  terminated  at  the  center  of  the  patch. 

It  should  be  noted  that  moving  the  center  of  the  patch  in  the  positive  y 
direction  results  in  a  different  impedance  than  when  it  is  moved  in  the 
negative  y  direction.  This  indicates  that  the  y  directed  components  of 
the  equivalent  magnetic  current  source  may  be  exciting  higher  order 
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nodes.  The  excitation  of  higher  order  nodes  would  cause  the  input  impedance 
to  be  different. 

It  is  desirable  fron  a  circuit  designer's  point  of  view  to  be  able 
to  predict  the  input  inpedance  of  the  antenna  so  that  the  appropriate 
matching  networks  can  be  designed.  This  is  especially  important  in 
integrated  "active"  array  applications.  The  cavity  model  is  a  good  means 
for  achieving  this  capability. 

Equivalent  magnetic  current  line  sources  can  be  placed  in  the  cavity. 
Modal  matching  techniques  can  then  be  used  to  find  the  resulting  fields 
in  the  cavity.  From  these  fields,  total  radiated  power,  dissipated  power 
and  stored  energy  can  be  calculated  in  order  to  predict  input  impedance. 

A  current  distribution  on  the  feed  needs  to  be  assumed.  The  current 
distribution  on  an  embedded  microstrip  feed  calculated  by  the  moment 
method  could  give  some  insight  into  the  electric  current  distribution  on 
the  center  conductor  of  the  coplanar  waveguide.  This,  in  turn,  could  be 
used  to  calculate  the  current  distribution  of  the  equivalent  magnetic 
current  source. 

At  this  point  a  computer  program  based  on  the  cavity  model  is  under 
development  to  carry  out  the  calculation  of  input  impedance. 

A  two-element  array  (Fig.  4)  was  constructed  and  tested  to  verify  the 
performance  of  the  coplanar  waveguide  feed  in  a  multielement  configuration. 
Radiation  pattern  plots  (Fig.  5  *  Fig.  6)  indicate  that  the  coplanar 
waveguide  feeds  are  coupling  power  into  the  antenna.  The  patterns  are 
also  symmetric  for  this  design.  More  work  needs  to  be  done  to  find 
the  optimum  feed  network  for  multiple-element  configurations. 
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Single  slot  lines  (Fig.  7)  can  also  be  used  to  feed  the  microstrip 


antenna.  To  excite  the  antenna,  it  is  necessary  to  orient  the  slot  in 
parallel  with  a  component  of  the  magnetic  field  of  the  mode  that  is  to  be 
excited.  Figures  S,  9  and  10  give  the  input  impedance  and  radiation  pattern 
data  for  a  4.1  GHz  slot  line  fed  antenna.  Again,  a  good  transition  from 
the  50  ohm  coaxial  cable  to  the  slot  line  was  constructed  so  that  the  measured 
input  impedance  of  the  slot  differs  from  the  actual  impedance  of  the 
antenna  by  only  a  phase  shift  factor.  A  1.1:1  input  VSWR  was  achieved  at 
resonance  indicating  that  virtually  all  the  incident  power  was  coupled  to 
the  slot  line  and  the  antenna. 

Several  observations  should  be  made  on  the  data.  The  H  plane  radiation 
pattern  is  asymmetric  and  the  impedance  locus  does  not  go  to  the  edge  of 
the  Smith  chart  when  the  frequency  is  moved  away  from  the  resonant  frequency. 
This  suggests  that  the  slot  line  itself  may  be  radiating.  Further  experiments 
will  be  required  to  verify  this. 

One  simple  experiment  that  was  performed  consisted  of  measuring  the 
return  loss  versus  frequency  for  both  the  slot  line  fed  antenna  and 
the  short-circuited  slot  line  (2  inches  long)  with  the  antenna  remov¬ 
ed  (Fig.  11).  For  the  short-circuited  slot  line  case,  a  return  loss  of 
approximately  3  dB  is  measured  at  4  GHz.  This  is  somewhat  greater  than 
what  one  would  expect  of  a  good  transmission  line  (i.e.,  low  loss)  with 
only  copper  and  dielectric  losses. 

One  possible  solution  to  the  radiation  problem  is  to  use  a  higher 
dielectric  constant  substrate,  which  will  tend  to  confine  the  fields. 

Another  possible  solution  is  to  configure  several  slot  line  fed  antennas 
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in  such  a  way  so  that  the  radiation  from  the  feeds  is  of  the  proper 
polarization  and  does  not  cause  an  asymmetry  in  the  radiation  patterns. 


CHAPTER  II. 

MEFSET  AMPLIFIER  DESIGN 


In  order  to  integrate  the  electroaagnetically  fed  microstrip  antenna 
with  a  MESFET  amplifier,  it  was  necessary  to  design  the  matching  network 
using  a  coplanar  transmission  line  rather  than  microstrip  line.  Because  of 
the  lack  of  research  done  in  the  literature  on  coplanar  discontinuities 
and  matching  networks,  there  were  neither  design  equations  (other  than 
equations  for  impedance  and  effective  dielectric  constant)  available  to  design 
appropriate  coplanar  elements  nor  computer-aided  design  packages  such  as 
SUPERCOMPACT  or  Touchstone  to  simulate  coplanar  structures.  Therefore, 
considering  the  current  lack  of  design  tools  as  well  as  the  narrow-band 
nature  of  the  patch  antenna,  a  simple  narrow-band,  4  GHz  maximum  available 
gain  (MAG)  amplifier  was  designed  rather  than  a  more  complicated  broadband 
amplifier  network  which  would  have  inevitably  required  some  computer-aided 
design.  In  addition,  the  transistor  was  assumed  to  be  unilateral  with 
source  and  load  impedances  equal  to  50  ohms  though,  in  general,  it  should 
be  noted  that  these  impedances  are  complex. 

The  design  of  the  MAG  amplifier,  which  was  integrated  with  the 
electromagnetically  fed  patch  (used  as  a  receiving  element),  can  be  summarized 
in  essentially  three  steps.  The  first  step  was  to  obtain  accurate  MESFET 
S-parameters  so  that  stability  and  conjugate  matching  impedances  can  be 
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determined.  The  second  step  was  to  determine  a  lumped  element  prototype 
which  achieves  the  conjugate  matching  (assuming  50  ohm  input  and  output 
impedances)  at  the  NESFET's  gate  and  drain.  The  third  and  final  step  was 
to  replace  the  lumped  elements  with  various  lengths  of  coplanar  transmission 
lines  to  approximate  the  lumped  elements. 

The  MESFET  used  in  the  amplifier  design  was  the  NEC  700  packaged 
transistor.  Though  manufacturer's  S-parameters  were  available,  it  was 
decided  to  measure  the  transistor's  S-parameters  due  to  the  sensitivity 
of  maximum  gain  amplifier  design  to  transistor  variation.  In  order  to 
measure  the  S-parameters  of  the  FET.  a  50  ohm  coplanar  thru  transmission 
line  section  was  fabricated  and  its  S-parameters  measured  on  the  network 
analyzer.  The  FET  was  then  placed  in  the  center  of  the  thru  line  and  then 
biased  at  a  drain  voltage  of  3  V  and  drain  current  of  30  mA  (corresponding 
to  the  Q  point  yielding  the  highest  gain).  The  S-parameters  of  FET/thru 
line  were  then  measured  on  the  network  analyzer.  The  S-parameters  of  the 
FET  were  calculated  from  the  FET/thru  line  data  by  deembedding  the  contri¬ 
bution  of  the  thru  line.  The  deembedding  was  accomplished  by  first  converting 
the  S-parameters  into  transfer  or  T-parameters  by  using  the  following 
relations : 
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Because  T-parameters  are  cascadable  under  matrix  multiplication,  we  were 
able  to  manipulate  the  T-matrlces  and  approximately  deembed  the  thru 
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line  fro*  the  FET/thru  line  measurements  and.  hence,  obtain  the  S-parameter 
characteristics  of  the  FET  only.  After  obtaining  the  desired  S-parameters , 
the  stability  of  the  FET  was  calculated  at  4  GHz  using  the  following 
stability  criteria 


K  = 


1  -|S„|2-  |Sj2|2+  |D|2 
2 1  S12| 


>  1  for  stablility 


where  D  =  |  S ^ |S  22"”  ^21^12! 


The  conjugate  impedances  were  also  calculated  from  the  Sll  and  S22  values 
of  the  FET  with  the  aid  of  a  Smith  chart  (16.5  «•  j66  ohms  and  52.5  - 
j82 . 5 ,  respectively). 

After  determining  the  stablity  and  the  required  conjugate  matching 
of  the  transistor  at  4  GHz.  a  lumped  element  matching  network  was  calculated 
using  standard  Smith  chart  techniques.  Though  there  are  a  multitude  of 
matching  schemes  one  could  use  at  this  stage,  two  of  the  most  basic 
matching  networks  were  developed  and  optimized  on  SUPERCOMPACT.  The  first 
design  (Fig.  12)  consisted  of  series  inductors  and  shunt  capacitors.  The 
second  design  (Fig.  13)  consisted  of  lengths  of  50  ohm  transmission  line 
and  series  inductors. 

Because  the  microstrip  patch  antenna  was  to  be  fed  electromagnetically 
via  coplanar  lines,  a  coplanar  distributed  matching  network  was  needed  in 
order  to  integrate  the  two  devices.  As  stated  earlier,  the  major  problem 
with  distributed  coplanar  networks  is  the  general  lack  of  design  equations 
with  which  to  approximate  lumped  elements.  There  are.  however,  semiempirical 
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design  equations  for  the  impedance  and  effective  dielectric  constant  of 
coplanar  lines.  Specifically,  these  equations  are 


£eff  *  K £r  +  ][  tanh  {  1.785  log  ( h/w)  +  1.75}  + 

{kW/h}  {.04  -  .7k  +  .01(1  -  .le0)(.25+k)}] 


impedance  of  line  = 


30  rc  K'(k) 
VEeff  K(k) 


where  K  =  elliptic  integral  of  the  first  kind 
K’(k)  =  K(k') 
k  =  — 

s  +  2w 

w  =  center  conductor  to  ground  spacing 
s  «  width  of  center  conductor 


The  only  matching  networks  that  could  be  constructed  with  such  limited 
design  equations  consist  of  high  impedance  and  low  impedance  coplanar 
transmission  lines.  It  is  well  known  that  a  length  of  high  impedance  line 
situated  between  two  relatively  low  impedance  lines  approximates  a  series 
Inductor  according  to  the  formula 


where  to  *  frequency  of  operation 
1  =  length  of  line 
z  =  impedance  of  line 
v  *  effective  velocity  of  line. 


Conversely,  a  low  impedance  line  situated  between  two  relatively  high 
impedance  lines  approximates  a  shunt  capacitor  according  to  the  formula 
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<0  C  = 


zv 


Though  it  was  theoretically  possible  to  realize  both  our  lumped  designs 
on  coplanar  transmission  lines,  it  turned  out  that  only  the  design 
consisting  of  the  series  inductors  and  lengths  of  transmission  line  was 
readily  realizable  in  coplanar  form.  This  result  was  due  to  the  fact  that 
low  impedance  lines,  i.e..  impedances  less  than  35  ohms,  which  are  needed 
to  approximate  shunt  capacitors,  are  very  difficult  for  us  to  fabricate 
especially  for  low  dielectric  substrates  where  the  center  conductor-to-ground 
distance  must  be  less  than  0.125  mm.  High  impedance  lines,  i.e.,  Impedances 
greater  than  100  ohms,  which  are  needed  to  approximate  series  inductors 
are,  on  the  other  hand,  easy  to  fabricate  due  to  the  relatively  large  center 
conductor-to-ground  distance.  Thus,  the  series  inductor/transmission  line 
network  was  chosen  to  be  realized  in  the  coplanar  structure.  The  choice 
of  the  high  impedance  value  to  be  used  was  somewhat  arbitrary.  Though  it 
is  known  that  the  inductor  approximation  improves  with  higher  impedance 
values,  we  note  that  an  upper  impedance  limit  exists  for  a  coplanar  line 
of  approximately  150  ohms.  That  is,  any  impedance  greater  than  150  ohms  will 
result  in  the  excitation  of  higher-order  non  TEM  modes.  With  this  limit 
in  mind,  a  high  impedance  value  of  120  ohms  was  selected.  After  having 
decided  what  impedances  to  use,  it  was  necessary  to  know  the  effective 
dielectric  constants  so  that  appropriate  physical  lengths  of  coplanar 
lines  could  be  calculated.  However,  because  the  design  equations  are 
empirically  corrected  for  coplanar  lines  on  high  dielectrics,  i.e.,  a  di¬ 
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electric  constant  greater  than  9,  rather  than  on  low  dielectrics,  which  is 
what  we  wished  to  use,  50  oha  and  120  oha  thru  lines  having  the  same  center 
conductor  width  of  3.5  ma  were  fabricated  on  62.5  ail  rexolite  (dielectric 
constant  of  2.62)  using  the  iapedance  design  equations.  They  were  then 
measured  in  the  laboratory  using  tiae  doaain  ref lectometry .  By  measuring 
the  unit  step  response  of  the  coplanar  lines,  we  were  able  to  determine 
the  reflection  coefficient  of  the  line  and.  hence,  determine  the  impedance 
of  the  line.  Then  using  the  following  relation, 

z0(  i  +r) 

( i  -  r) 

where  Z0=  50  Cl  =  characteristic  impedance  of  source 
r  *  reflection  coefficient , 

we  can  deteraine  the  iapedance  of  the  line.  According  to  the  step  response 
(Fig.  14)  the  reflection  coefficient  of  the  50  oha  coplanar  line  equals 
-0.02  which  corresponds  to  an  iapedance  of  48  ohas.  The  reflection  coefficient 
of  the  120  oha  coplanar  line  (Fig.  15)  was  found  to  be  0.4  which  corresponds 
to  an  impedance  of  117  ohas.  Therefore,  the  iapedance  design  equation  appeared 
to  be  valid  for  low  dielectric  substrates.  To  measure  the  effective 
dielectric  constant  of  the  lines,  we  needed  to  aeasure  the  impulse  response 
of  the  lines  with  short  circuits  placed  at  fixed  distances  along  the  line.  By 
knowing  the  impulse  response  of  the  line,  we  could  aeasure,  in  tiae,  how 
far  apart  the  short  circuits  were.  Since  the  distances  between  the  short 
circuits  are  known  (2  cm),  we  were  able  to  calculate  the  velocity  of 
the  impulse  and,  hence,  the  effective  dielectric  constant  of  the  line.  A- 
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nalysis  of  TDR  measurements  shows  that  the  effective  dielectric  constants 
were  1.485  and  1.395  for  the  50  oha  and  120  ohm  lines,  respectively.  Comparing 
the  calculated  dielectric  constants  of  1.695  for  the  50  oha  case  and 
1.421  for  the  120  oha  case  with  the  TDR  aeasureaents ,  we  see  good  agreement 
for  the  120  ohm  case  but  a  rather  large  discrepancy  between  the  effective 
dielectric  constants  for  the  50  ohm  case.  The  fact  that  the  measured 
effective  dielectric  constants  were,  in  general,  less  than  those  predicted 
by  the  design  equations  indicates  that  the  field  confinement  for  the  low 
dielectric  case,  particularly  for  the  narrow  center  conductor-to-ground 
spacing,  is  less  than  that  predicted  by  the  design  equations.  It  is 
believed  that  this  difference  is  caused  primarily  by  the  empirical  corrections 
added  to  the  effective  dielectric  equations.  Since  the  equations  were 
made  to  fit  experimental  data  obtained  for  coplanar  lines  on  high  dielectric 
substrates  rather  on  than  the  low  dielectric  substrates,  it  is  a  likely 
source  of  error.  For  our  coplanar  design,  we  used  the  experimental  effective 
dielectric  constants  to  convert  the  transmission  line/series  inductor  MAG 
amplifier  matching  network  into  lengths  of  50  ohm  and  120  ohm  coplanar  lines. 

CHAPTER  III. 

INTEGRATED  "ACTIVE"  ARRAY  ELEMENT 

After  designing  the  transmission  line/series  inductor  amplifier  on 
coplanar  transmission  lines,  a  prototype  was  constructed  on  62.5  mil 
Rexolite  (Fig.  18).  The  amplifier  was  then  tested  on  the  network  analyzer 
and  found  to  have  a  return  loss  of  16.7  dB  and  gain  of  9  dB  at  a  frequency 
of  3.78  GHz  (Fig.  16).  Though  the  frequency  of  operation  is  slightly  off 
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the  design  frequency  of  4  GHz  and  the  gain  is  approximately  5  dB  less 
than  the  maximum  available  gain,  the  coplanar  matching  network  has  nonetheless 
improved  the  performance  of  the  PET.  The  deviation  in  performance  from 
theoretical  performance  is  due  to  the  presence  of  coplanar  discontinuities, 
dc  blocking  capacitors  as  well  as  the  biasing  network  (fig.  18),  all  of 
which  were  either  currently  unmodelable  on  the  computer  or  assumed  to 
have  negligible  effect  on  the  matching  network. 

After  measuring  the  amplifier  by  itself,  the  amplifier  was  integrated 
with  a  single  microstrip  patch  (operating  frequency  of  3.78  GHz).  The 
patch  was  placed  on  the  uncladded  side  of  the  Rexolite  substrate  over  a 
length  of  50  ohm  coplanar  line  which  extended  from  the  matching  network 
(fig.  18).  In  order  to  satisfy  the  50  ohm  source  assumption  used  in  the 
amplifier  design,  the  input  impedance  of  the  patch  antenna  at  various 
coplanar  feed  points  was  determined  by  measuring  the  return  loss  of  the 
patch  versus  feed  location  (Table  1).  A  feed  location  with  an  optimum 
return  loss  of  16.0  dB  was  found.  Since  such  a  return  loss  did  not  represent 
a  perfect  50  ohm  line,  a  mismatch  loss  given  by 


Mismatch  Loss  =  10  log 


(i-lrd2)d-lrsl2) 

H  -I  £11  Is  I  exPJ(  0L+  0s  )l2 


] 


where  |  j^|  *  reflection  coefficient  at  antenna 

|  jg  |  =  relection  coefficeint  of  the  amplifier 

Pdei  =  power  delivered 
Pav  =  power  available 
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was  expected  to  reduce  the  gain  of  the  amplifier. 


The  radiation  pattern  for  the  integrated  "active"  element  is  given 
in  Pigs.  19.  20.  It  is  superimposed  on  the  radiation  pattern  of  a  "passive" 
element  with  the  same  feed  structure  but  no  amplifier.  The  difference  in 
power  levels  indicates  that  the  amplifier  adds  approximately  13  dB  of 
gain  of  the  antenna  element.  The  ripple  in  the  radiation  pattern  of  the 
integrated  element  is  due  to  the  effects  of  diffraction  at  the  edges  which 
were  magnified  by  the  use  of  a  large  but  finite  ground  plane  on  which  the 
the  element  was  fabricated.  A  smaller  ground  plane,  which  reduces  the 
diffraction  effects,  was  used  for  the  pattern  measurement  of  the  antenna 
element  alone  and  found  to  reduce  the  ripple.  One  notes  that,  unlike  the 
H-plane  pattern,  the  E-plane  radiation  pattern  demonstrates  a  slight 
asymmetry.  This  asymmetry  is  attributed  to  the  presence  of  the  coplanar 
feed  and  amplifier  which  lie  on  the  same  axis. 


In  the  future,  work  will  be  focused  on  refining  the  integration 
design  procedure  by  eliminating  the  need  to  resort  to  empirical  design 
data  while  improving  the  overall  performance  characteristics  of  the 
integrated  antenna/amplifier.  In  particular,  a  computer  program  based  on 
theoretical  analysis  will  be  developed  which  will  calculate  the  input 
impedance  of  the  coplanar  and  slot  line  feed  patch  antennas.  Work  on 
characterizing  novel  coplanar  line  configurations  will  also  be  done  so 
that  a  more  compact,  broader  band  amplifier  can  be  constructed  which  will 
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not  only  ainiaize  the  asyaaetry  in  the  radiation  pattern,  but  facilitate 
aultieleaent  antenna/aaplifier  integration.  Finally,  attention  will  be 
focused  on  designing  integrnble  coplanar  feed  networks  which  can  be 
used  in  the  fabrication  of  large  aultieleaent  antenna/amplifier  arrays. 
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TABLE  1 


INPUT  REFLECTION  COEFFICIENT  VERSUS  FEED  LOCATION  FOR  COPLANAR  FEED 

MICROSTRIP  ANTENNA 


Ax 

(cm) 

Ay 

(cm) 

201og(I*.| 

(db) 

0 

1.000 

■4.S 

0 

0.875 

-4.9 

0 

0.750 

-7.9 

0 

0.500 

-11.6 

0 

0 

-16.0 

0 

-.0250 

-11.8 

0 

-0.500 

-5.8 

0 

-1.000 

-1.6 

0.250 

0 

-13.6 

0.500 

0 

-142 
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Figure  1.  Microatrip  antenna  eleoent  with  coplanar  waveguide  feed. 

Electric  Walls: 
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Figure  2.  Cavity  aodel  representation  of  nicroatrip  antenna  element  with 
coplanar  waveguide  feed. 
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180°  90°  0°  90°  180° 

ANGLE  (DEG) 

Figure  5.  Two  element  array  radiation  pattern  -  E-plane.  (Coplanar  waveguide 
feed)  (f  -  4.5  GHz) 


180°  90®  0  90®  180® 

ANGLE  (DEG) 

Figur#  6.  Two  <leunc  array  radiation  pattarn  -  H-plane.  (Coplanar  waveguide 
feed)  (f  •  4.5  GHz) 
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Figure  8.  Microstrip  antenna  eleaent  with  slot  line  feed. 

a)  Element  input  impedance  (z)  versus  frequency 

b)  Element  input  VSWR  versus  frequency. 


Figure  9.  Radiation  pattern  of  element  with  slot  line  feed.  (E-plane) 
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Figure  10.  Radiation  pattern  of  element  with  slot  line  feed.  (H-plane) 
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Figure 
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50 Cl  Load 


Figure  12.  Amplifier  design  #1. 


Figure  13.  Amplifier  design  #2. 
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Figure  16.  Frequency  response  of  prototype  aoplifier  #1. 
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Figure  20.  Radiation  patterns  for  "active”  and  "passive”  microstrip  antenna 
elements.  (E-plane) 
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CHAPTER  IV 


RIGOROUS  ANALYSIS  OF  A  CIRCULAR  PATCH  ANTENNA 
EXCITED  BY  A  MICROSTRIP  TRANSMISSION  LINE 

M.  Davidovitz 

Y.T.  Lo 
ABSTRACT 

Circular  patch  antenna  excited  by  a  microstrip  transmission  line  is  analyzed.  Boundary  con¬ 
ditions  are  enforced  on  a  portion  of  the  microstrip  feed  line  as  well  as  the  patch  antenna.  The 
integral  equation  for  the  unknown  currents  on  the  antenna  and  the  feed  is  solved  by  applying  the 
Galerkin  method  of  moments  in  the  Fourier  Transform  domain.  A  method  designed  to  accelerate 
the  convergence  of  the  resulting  spectral  integrals  is  outlined  and  discussed.  Extensive  input 
impedance  data  are  presented  for  the  case  of  the  eiectromagnetically  coupled  circular  patch.  The 
accuracy  of  the  computed  results  is  verified  experimentally. 
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I.  INTRODUCTION 


Microstrip  antennas  have  been  analyzed  extensively  in  the  past  decade,  employing  a  wide 
variety  of  analytical  techniques,  ranging  from  the  simple  intuitive  models  to  rigorous  numerical 
treatments.  Oftentimes,  in  order  to  make  the  computations  more  tractable,  simplifying  assump¬ 
tions  about  the  feed  were  introduced  into  the  formulation.  This  placed  restrictions  on  the  range  of 
validity  of  the  models  in  question.  Specifically,  probe  and  microstrip  transmission  line  feeds  were 
in  many  instances  modeled  by  current  filaments. 

Aside  from  the  fact  that  such  idealizations  are  not  adequate  for  thicker  elements,  they  cannot 
be  applied  for  more  complicated  feeding  arrangements.  Particularly,  in  the  case  of  the  electromag¬ 
netic  coupling  or  proximity  feed  it  is  clear  that  boundary  conditions  on  at  least  a  portion  of  the 
microstrip  line  must  be  accounted  for  in  the  formulation. 

This  paper  undertakes  the  analysis  of  the  circular  patch  antenna  excited  by  a  microstrip 
transmission  line.  The  method  presented  is  applicable  to  all  microstrip  line  feeding  arrangements, 
including  the  edge  feed,  capacitive  gap  feed  and  the  proximity  or  coupling  feed.  The  problem  is 
formulated  in  the  Fourier  Transform  (FT)  domain.  Moment  Method  is  applied  to  solve  the  integral 
equation  for  the  currents  on  the  patch  and  portion  of  the  feed  line.  Considerable  attention  is  given 
to  acceleration  of  convergence  of  the  Sommerfeld-lype  integrals  which  arise  in  the  course  of  the 
analysis. 

In  Section  II  the  integral  equation  for  the  unknown  patch  and  microstrip  feed  line  currents  is 
formulated.  The  derivation  is  performed  in  the  FT  domain  and  utilizes  the  dyadic  Green  function 
for  the  multilayered  grounded  dielectric  slab.  A  set  of  basis-testing  functions  is  introduced  for  the 
patch  and  microstrip  line.  Parseval's  theorem  is  applied  to  implement  the  Galerkin  Moment 
Method  solution  in  the  FT  domain.  The  resulting  matrix  elements  take  the  form  of  improper  single 
or  double  integrals. 

Section  III  outlines  and  discusses  a  method  to  accelerate  the  convergence  of  the  improper 
integrals.  The  idea  behind  the  method  is  similar  to  singularity  subtraction  in  evaluation  of  singular 
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integrals.  It  involves  subtraction  of  the  slowly  decaying  terms  of  the  integrand,  which  due  to  their 
simple  form  can  be  readily  evaluated.  This  procedure  tends  to  increase  the  efficiency  and  accuracy1 
of  the  computation. 

The  degree  of  accuracy  to  which  the  input  impedance  can  be  computed  depends  to  a  large 
extent  on  the  model  of  the  excitation  mechanisms.  In  Section  IV  a  unit  voltage  source  is  introduced 
on  the  microstrip  transmission  line.  Input  impedance  is  defined  based  on  the  transmission  line 
model  of  the  actual  physical  configuration. 

Theoretically  computed  results  are  compared  with  experimental  data  in  Section  V.  Since  vir¬ 
tually  no  data  exists  on  proximity  fed  circular  patches,  efforts  are  concentrated  in  that  direction. 
Results  are  exhibited  for  a  number  of  configurations. 

Finally,  discussion  and  conclusions  are  presented  in  Section  VI. 


1  Note,  the  ultimate  test  of  accuracy  is  comparison  of  computed  results  with  experimental  data. 
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n.  FORMULATION  OF  THE  PROBLEM 


There  exist  several  equivalent  ways  of  deriving  the  dyadic  Green  function  for  the  multilay¬ 
ered  grounded  dielectric  slab.  The  present  approach  utilizes  the  modal  decomposition  of  tie  fields 
transverse  to  z  into  E(T\1)  and  H(TE)  components.  The  original  three-dimensional  problem  is 
thereby  reduced  to  a  one-dimensional  transmission  line  problem.  This  procedure  is  detailed  else¬ 
where  [l ]:  the  resulting  equations  are  summarized  below. 

To  facilitate  the  following  presentation  of  the  dyadic  Green  function,  several  necesarv 
definitions  are  introduced.  The  FT  relations  used  throughout  the  paper  are: 

uo 

F{A(p.z)|  =  A(k,.z)  =  JJ  A(p.2)ejk|P  dp  (la) 


oo 

F"1lA(k,.z)}  =  A(p.2)  =  ■  ff  A(kt.z)e_ik,p  dk, 

v  ZTT  )  — oo 


(lb) 


where  k,  is  the  transverse  vector  wavenumber,  e.g..  in  the  Cartesian  coordinate  system 

k,  =  kxx  +  kyy  =  k,(cosa  x  +  sina  y)  ;  k,2  =  kx2  +  ky2  (2) 

The  corresponding  unit  vector  will  be  denoted  by  k,. 

Consider  a  transverse  electric  current  distribution  situated  in  a  stratified  medium  with  layers 
perpendicular  to  the  z  axis.  The  FT  of  the  dyadic  Green  function  for  this  configuration  takes  on  the 
following  form: 


G(k,;z.z')  -  — V'(z.z')k,k,  -  V''(z.z')(kIxz)(klxz)  +  —  I'(z.z')zk, 

dt€ 


(3) 


V'(z.z  )  and  V”(z.z  )  are.  respectively,  the  E-  and  H-mode  characteristic  Green  functions, 
satisfying  the  following  set  of  transmission  line  equations 


V(z.z')  +  jk2ZI(z.z')  =  0 
dz 


(4a) 


— I(z.z')  +  jk2YV(z.z')  = 
dz 


(4b) 
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where  i  -  5(z— z)  or  i-0  depending  on  whether  the  source  point  z  and  the  observation  pont  z  are  or 
are  not  in  the  same  stratum.  Furthermore: 

Z‘  =  1/Y’  =  kz/<ve  ,  E— modes  (5a) 

Z  = 

Z"  =  1/Y"  =  (Dfi/  kz  .  H— modes  (5b) 

kz  =  *v/k2— kt2:  k2  =  (*>2a ie  (5c) 

The  branch  of  the  square  root  in  equation  (5c)  is  chosen  such  as  to  satisfy  the  radiation  condition 
for  the  assumed  e^1  time  dependence. 

As  an  example  illustrating  the  solution  of  equations  (5a.b)  consider  the  case  of  a  two-layer 
slab  backed  by  a  ground  plane.  Figure  1  shows  the  physical  configuration  alongside  the  equivalent 
transmission  line  model.  The  dielectric  layers  are  characterized  by  their  permittivities  6j.  €2  and 
respective  thicknesses  tj.  t2.  The  source  point  z  is  located  in  the  second  layer,  i.e.  z^z'^zj.  The 
characteristic  Green  function  for  this  example  is  given  by: 


V,(z2.z')  e-^1-^ 

z2  <  z 

(6a) 

V(z.z')  =  V2(z.z') 

z,^z^z2 

(6b) 

sin(kllz) 

V  2(Z,.Z  )  — : — 7- - r 

sin(kjiZi) 

O^Z^Z1 

(6c) 

where 

V2(z.z’)  =<=5-^ - [coskz2(z<— zj  +  jY'(z|)sinkz-,(z<— Zj)]  x 

Y  (z,) 

X  (coskz2(z>-z1)— f?'(z,)sinkl2(z>-z1)]  (7a) 
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\_  ^2—j^nCOl(kj2l2) 
Ul>=  Y(1-jY2cot(kz2i2) 

Y(z,)  *  “j-rA  COt(kj,t,) 
>2 


:  ^f(zj)  =  Y2^'(Z|) 

^(2^  =  y2V(2,) 


Y  (z,)  s^fz,) +"?(z,). 


(7b) 


(7c) 


and 


z,  z  <  z' 

Z*=  .  .> 
z  .  z  <  z 

Note,  equations  (4.6.7)  are  equally  valid  for  both  E-  and  H-mode  characteristic  Green  func¬ 
tions.  V'(z.z  )  or  V'  (z.z  )  are  obtained  from  (6.7)  by  substitution  of  the  respective  modal  admit¬ 
tances  Y  or  Y”.  given  by  (5). 


Formulation  of  the  integral  equation  for  the  unknown  currents  on  the  patch  and  feedline.  J,, 
and  J,.  respectively,  can  be  carried  out  with  the  aid  of  the  Green  function  given  by  equations  (3-7). 
Enforcement  of  the  boundary  condition  requiring  the  total  tangential  electric  field  to  vanish  on  the 
microstrip  patch  and  line,  yields  the  sought  integral  equation 


z  X 


E'(r)  +  E(r)  =0  re  5.  +  S, 


where  E'(r)  denotes  the  impressed  or  incident  electric  field  and 

E(r)  -  f  f  G(r x)  J,(r')dr'  +  f  f  G(r.r’)  J,(r')dr' 

s,  S, 


(8a) 


G(k(;z.z')  •  Jf,(kt) 


+  F~l 


G(k,:z.z')  •  J,(kt) 


(8b) 


Equation  (8b)  is  obtained  from  (8a)  by  extending  the  domains  of  integration  Sp  and  S,.  shown 
in  Figure  2.  to  the  entire  plane  and  applying  basic  FT  relations. 


Galerkin  s  method  of  moments  is  implemented  in  the  FT  domain  to  reduce  the  integral  equa¬ 
tion  (8b)  to  a  matrix  equation.  The  accuracy  and  efficiency  of  the  solution  largely  depends  upon 

i 

i 
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the  choice  of  the  expansion  and  testing  functions.  This  choice  can  be  guided  by  previous  investiga¬ 
tions  of  related  geometries,  and  dictated  by  the  inherent  features  of  the  configuration,  such  as  sym¬ 
metries.  edges,  etc.,  as  well  as  numerical  considerations.  The  selection  of  the  expansion-testing 
functions  for  the  problem  at  hand  is  based  upon  the  previous  work  of  several  authors  [2.3.4], 

The  currents  on  the  circular  patch  of  radius  c  are  represented  by  the  following  set  of  basis 
functions: 


J/>  =  Z  lm*  J<m-  +  Z  14  m.n.n.£  =  0.1.2 .  m  +  v.n  +  |  odd  (9*) 

m  v  n.£ 

where  J^,,  =  pP^p/ckosfm#):  P,,(p/c)  =  L'v(p/c)"v/ 1— (p/c)2  (9b) 

<t>£(p/c)sin(n<£):4>£(p/c)  =  Tf(p/c)/V  1—  (p/c.)2  (9c) 

T^.  U„  are  Cheybishev  polynomials  of  first  and  second  kind,  respectively.  These  basis  functions  are 
particularly  suitable  for  problems  involving  the  circular  disc  because  they  correctly  model  the 
behavior  of  the  currents  near  its  edge.  In  addition  they  possess  a  closed  form  Fourier-Bessel 
transforms,  which  are  given  by 


Jp 

•n u/ 


=  7rjm+1 


k,  cos( ma)U~v( k tc/2 )  —  (kt  x  z)sin(ma  )£/,£*(  k,c/2) 


(10a) 


Jnl  =  "T+l 


kt  cos( no )T*fi k vc/2 )  —  (lct  x  z)sin(na)7’^(k,c/2) 


(10b) 


Functions  (/„„( k,c/2).  T^( k,c/2)  can  be  found  in  Appendix  A. 

The  current  on  the  microstrip  transmission  line  of  width  2w  is  approximated  by  a  subsec¬ 
tional  basis  set  consisting  of  the  following  functions 

l  =  xZ  W  =  xZ  lnT„(x)P(2w)  ( 1  la) 

n  n 


where 
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Tn(x)  = 


I  x-xn  I 
A 


xn— A  ^x^xn+A 


(Ub) 


P(2w)  = 


2^7  — w  <  y  ^  w 

0  otherwise 


with  w  as  shown  in  Figure  2.  The  FT  of  Jnx  is  readily  shown  to  be 


(11c) 


Ax  =  A  « 


J*X*«  — 


sin 4 


kxA 


k.A 


sin(kyw) 

(kyw) 


(12) 


Next,  the  moment  method  matrix  is  formulated  using  the  basis-testing  functions  introduced 
above.  Both  expansion  and  testing  take  place  in  the  FT  domain,  the  latter  with  the  aid  of  the  fol¬ 
lowing  Parseval's  Theorem 


co  , « 

fj  E(p.z)  -J(p.z)ds  =  <»> 

The  resulting  matrix  can  be  partitioned  into  three  distinct  submatrices.  One  of  the  subma¬ 
trices.  itself  of  a  block  diagonal  form,  represents  the  interaction  between  the  various  current  modes 
on  the  patch.  Another  serves  the  same  purpose  for  the  microstrip  line  and  possesses  a  Toeplitz 
structure.  The  third  couples  the  patch  current  modes  to  the  microstrip  line.  Tabulation  of  the 
various  types  of  elements  can  be  found  in  Table  I. 

To  recapitulate,  the  Green  function  for  the  stratified  dielectric  slab  backed  by  a  ground  plane 
was  used  to  formulate  the  integral  equation  for  the  currents  on  the  patch  antenna  and  microstrip 
feedline.  A  set  of  expansion-testing  functions  was  introduced  and  Galerkin  s  method  of  moments 
applied  to  reduce  the  problem  to  a  set  of  linear  equations. 
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IE.  NUMERICAL  ANALYSIS 


As  Table  I  reveals,  the  moment  method  matrix  elements  are  given  by  improper  single  and 
double  integrals  of  very  complex  integrands.  Evaluation  of  such  integrals  requires  brute  force 
numerical  integration  and  as  such  tends  to  be  very  expensive.  The  accuracy  and  efficiency  of 
integration  can  be  significantly  improved  by  a  procedure  presented  in  this  section.  The  method  is 
analogous  to  the  singularity  subtraction  method  in  evaluation  of  singular  integrals  and  to 
Rummer  s  transformation  of  slowly  convergent  series. 


Implementation  of  the  accelerating  procedure  requires  the  prior  knowledge  of  the  asymptotic 
expansions  of  V".  V"  as  k,  —  oo.  Table  II  gives  the  first  term  of  the  asymptotic  expansions  of 
V'\  V”  for  a  two  layer  grounded  dielectric  slab.  All  the  relevant  dimensions  are  indicated  in 
Fig.  1. 

A  glance  at  Table  II  reveals  that  for  the  case  where  z  5*  z.  the  characteristic  Green  functions 
decay  exponentially  for  large  values  of  kt.  This  means  that  for  electromagnetically  coupled 
patches,  where  the  antenna  and  the  feed  are  vertically  separated,  the  integrands  of  Z£l  p.  Z^|p 
become  vanishingly  small  outside  of  some  finite  interval.  This  interval  rapidly  decreases  as  the 
antenna-feed  separation  is  increased.  However,  in  cases  where  the  source  and  the  observation  plane 
coincide,  i.e.  z  »  z\  the  integrands  exhibit  a  much  slower  algebraic  decay.  In  those  situations  it  is 
most  expedient  to  implement  the  accelerating  procedure  in  order  to  increase  the  efficiency  of 
integration. 

A  brief  illustration  of  the  accelerating  technique  now  follows.  Consider  one  of  the  matrix 
elements,  namely  Z£.  found  in  Table  I.  Letting  VA'  denote  the  first  term  in  the  asymptotic  expan¬ 
sion  of  V'.  the  expression  for  is  rewritten  as  follows: 


oo  w/2 


-pq 


--A// 


o  0 


[V‘  —  VA']  cosher  +  V  "sin2ar 


/2cos(k,  I  xp— xq  I  coso)k,dktda  + 
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+  — —f  f  V^cos2®  /x:cos(k,  I  xp— xq  I  cosa)k,dk,da 
V2  O  I) 

=Z“+ A«  (14) 

Note,  the  most  dominant  asymptotic  term  of  the  original  integrand  has  been  subtracted. 
Therefore,  for  large  values  of  kt.  the  modified  integrand  will  decay  considerably  faster.  The  inter¬ 
val  of  integration  in  k,  required  to  achieve  relative  convergence  is  thereby  reduced  from  —  150  k0  - 
200  k0  to  —  40  k0  -  50  k„,  where  k0  is  the  free  space  propagation  constant.  The  integral  denoted  by 
AjJ,*  has  a  much  simpler  form  and  can  be  evaluated  very  efficiently.  In  many  instances,  including 
the  case  considered  above,  the  asymptotic  integral  can  be  evaluated  in  closed  form  in  terms  of  ele¬ 
mentary  functions. 

Additional  numerical  analysis  can  be  performed  on  the  integrals  listed  in  Table  I.  For  exam¬ 
ple.  more  than  one  asymptotic  term  can  be  subtracted  from  the  original  integrand.  However,  the 
additional  asymptotic  terms  give  rise  to  singular  behavior  at  k,  *  0.  This  difficulty  can  be  circum¬ 
vented  by  introducing  a  positive  regularizing  parameter  u  and  replacing  k,  by  Vu2  +  k,2  in  the 
asymptotic  expansions.  This  technique  was  used  to  accelerate  slowly  convergent  infinite  series  in 
the  study  of  infinite  periodic  structures  [4.5]. 

The  integrands  can  be  further  simplified  by  extracting  the  surface  wave  pole  contributions 
and  integrating  them  analytically.  However,  since  in  all  practical  situations  the  dielectrics  used  are 
lossy,  the  surface  wave  poles  reside  off  the  real  axis  path  of  integration  and  present  no  major 
difficulty  in  numerical  integration. 

To  summarize,  a  method  accelerating  the  convergence  of  the  improper  spectral  integrals  was 
outlined  and  illustrated.  Application  of  this  method  increases  the  efficiency  of  matrix  elements 
evaluation.  Detailed  analysis  of  all  matrix  elements  using  the  accelerating  technique  can  be  found 
in  Appendix  B. 
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IV.  EXCITATION  MECHANISM  AND  INPUT  IMPEDANCE 


As  is  pointed  out  in  reference  [7],  one  of  the  most  difficult  aspects  of  the  problem  under  con¬ 
sideration  is  the  lack  of  a  suitable  model  for  the  excitation  mechanism.  In  most  practical  situations 
there  exists  a  transition  between  the  coaxial  cable  connector  and  the  microstrip  line  feeding  the 
antenna.  To  determine  the  true  impedance,  as  would  be  measured  by  the  network  analyzer,  the 
analytical  model  should  in  principle  account  for  the  coax  to  microstrip  transition  and  use  the  coax¬ 
ial  line  1  EM  mode  as  the  primary  excitation.  However,  this  kind  of  rigorous  analysis,  though  pos¬ 
sible  in  prinicple.  would  be  prohibitively  complicated  and  expensive  in  computation. 

References  [7.8.9]  advance  several  excitation  models  based  upon  the  assumption  that 
sufficiently  far  from  the  antenna  the  microstrip  line  supports  only  the  quasi-TEM  mode.  A 
transmission  line  analogy  of  the  feed  is  then  introduced,  wherein  the  antenna  is  seen  as  an 
impedance  termination.  The  transmission  line  is  excited  by  a  lumped  voltage  source,  as  shown  in 
Figure  3.  The  characteristic  impedance  Zc  and  the  propagation  constant  0  can  be  calculated  using 
the  standard  formulas  for  microstrip  lines,  although  the  presence  of  a  dielectric  overlay  will 
modify  these  quantities  somewhat. 

In  the  context  of  the  moment  method  solution,  the  lumped  voltage  excitation  of  the  transmis¬ 
sion  line  is  equivalent  to  the  della  gap  generator.  Therefore,  the  excitation  vector  [Vj  has  the  fol¬ 
lowing  simple  form 


V,  = 


1  .  at  the  position  of  the  generator  (x,  =  x  ) 
0  .  otherwise 


(15) 


One  of  the  major  advantages  of  this  excitation  model  is  clearly  the  simplicity  of  the  voltage  vector 
evaluation. 

Solution  of  the  matrix  equation  for  the  currents  on  the  patch  and  the  feed  line  yields  the 
value  of  the  current  I(x  )  at  the  location  of  the  generator.  A  straightforward  analysis  of  the 
configuration  in  Figure  3  shows  that 
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(16) 


=*Z(x')  +l(x) 

As  Fig.  3  indicates,  the  left  end  of  the  transmission  line  is  terminated  in  an  open  circiafc.  This 
fact  is  a  consequence  of  chocsing  basis  functions  which  force  the  current  to  zero  at  the  end  of  the 
microstrip  line.  Hence  Z(x')  is  equal  to  the  impedance  of  a  section  of  open  circuited  transmission 
line  of  length  (x^  —  x)+A where  A accounts  for  the  end  capacitance  [10],  The  antenna  input 
impedance^(x')  can  now  be  computed  using  (16)  in  the  following  manner: 

ZIN  =l(x')  =  -’Z(x') 

where  v(x')  =  —1.  and  I(x’).  and*Z(x‘)  are  known  quantities. 
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V.  THEORETICAL  AND  EXPERIMENTAL  RESULTS 


The  theory  outlined  thus  far  is  sufficiently  general  to  allow,  by  slight  adjustments  of  parame¬ 
ters.  analysis  of  most  microstrip  transmission  line  feeding  arrangements.  Those  include  the  direct 
edge  connection,  capacitive  gap  feed  and  the  two-laver  or  proximity  feed.  The  latter  case  is  partic¬ 
ularly  interesting  [7.11-13]  because  it  allows  the  transmission  line  to  be  placed  close  to  the  ground 
plane,  thereby  reducing  undesirable  radiation  from  the  feed  network.  At  the  same  time  the 
bandwidth  of  the  microstrip  antenna  can  be  enhanced  by  increasing  the  patch  to  ground  spacing. 
The  results  presented  in  this  section  will  also  show  that  the  proximity  feed  has  an  additional 
advantage  over  other  microstrip  feeds  in  that  a  perfect  match  to  a  50  SI  line  can  be  achieved  simply 
by  selecting  the  proper  coupling  distance  d  (see  Fig.  2). 

In  this  section  computed  input  impedance  results  for  the  case  of  the  electromagnetically  cou¬ 
pled  circular  patch  are  compared  with  experimental  data.  Shown  in  Figs.  4.5  are  input  impedance 
loci  for  a  patch  of  radius  17.5  mm.  The  dielectric  portions  of  the  antenna  are  composed  of  Rexolite 
2200  with  relative  permittivity  of  2.62  and  a  dielectric  loss  tangent  of  approximately  .001.  In  all 
cases  the  microstrip  line  exciting  the  antenna  is  4.6  mm  wide  and  is  placed  1.6  mm  above  the 
ground  plane.  The  impedance  reference  plane  is  in  all  cases  located  80  mm  away  from  the  center  of 
the  patch. 

The  computed  results  in  Figs.  4.5  were  obtained  utilizing  the  moment  method  approach  out¬ 
lined  in  Section  II.  Currents  on  the  patch  were  approximated  by  entire  domain  basis  functions, 
given  by  (9).  Compared  with  a  subsectional  expansion,  fewer  basis  functions  were  required  to 
obtain  numerical  convergence,  resulting  in  a  more  efficient  solution.  The  results  in  Figs.  4.5  were 
computed  using  the  basis  functions  in  Eqs.  (9a)-(9c)  with  the  following  index  set:  {[m.n-0;  v,£  «= 
1.3].  [m.n-I;  m  0.2].  [m.n=2;  *  1.3.5],  [m.n  =  3:  v.£  =  0.2. 4.6]).  The  currents  on  the  micros¬ 

trip  transmission  line  were  expanded  in  a  subsectional  basis  set  and  approximately  12-15  basis 
functions  per  guide  wavelength  were  employed  in  computing  the  data. 
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Figs.  4.5  display  good  agreement  between  theory  and  experiment.  The  discrepancies  which  do 
exist  can  be  attributed  to  both  theoretical  and  experimental  error.  The  inaccuracies  in  computed 
data  are  possibly  due  to  the  simplifying  assumptions  made  in  the  course  of  the  formulation.  For 
instance,  the  transverse  current  variation  on  the  microstrip  line  was  assumed  to  be  uniform  (see 
eqn.  (11)).  disregarding  the  singular  edge  behavior.  Based  on  the  analysis  of  the  transverse  current 
variation  on  printed  dipoles  performed  in  [5],  use  of  a  more  elaborate  transverse  variation  was 
deemed  unnecessary  for  line  width  considered  here.  One  of  the  sources  of  experimental  error  was 
the  discontinuity  encountered  in  the  transition  between  the  coaxial  connector  and  the  microstrip 
line. 


VL  DISCUSSION  AND  CONCLUSIONS 

A  method  for  computing  the  input  impedance  of  a  circular  patch  antenna  excited  by  a  micros¬ 
trip  transmission  line  has  been  presented  in  this  paper.  A  distinguishing  feature  of  the  present 
analysis  is  the  incorporation  into  the  formulation  of  the  boundary  conditions  on  the  feed.  This- 
type  of  a  rigorous  treatment  has  several  advantages,  namely  it  yields  accurate  input  impedance 
results,  and  is  applicable  to  a  wide  variety  of  feeding  arrangements  utilizing  the  microstrip 
transmission  lire.  The  negative  aspect  of  the  rigorous  treatment  is  the  additional  computational 
expense  incurred  as  a  result  of  including  the  feed  in  the  formulation.  Consequently,  the  efficiency 
of  numerical  computation  is  a  very  important  consideration  in  the  overall  solution.  This  aspect  of 
the  problem  was  addressed  in  Section  III.  where  a  method  of  increasing  the  efficiency  and  accuracy 
of  numerical  computation  was  presented. 

A  difficulty  of  a  theoretical  nature  encountered  in  the  analysis  was  the  absence  of  a  general 
excitation  model.  In  this  paper  the  input  impedance  was  defined  and  computed  using  an  excitation 
model  based  upon  the  assumption  that  the  microstrip  line  supports  a  quasi-TEM  mode.  This 
assumption  seems  to  hold  best  for  microstrip  lines  which  have  a  relatively  large  width  to  height 
above  ground  ratios.  Possible  recommendations  for  future  research  include  development  of  a  more 
general  excitation  model  and  analysis  of  more  sophisticated  coupling  structures. 
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APPENDIX  A 


The  functions  V„ 7t(k,c/2).7",*f  (k,c/2)  are  given  by 
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where  in  (A1-A4)  "J*  denotes  the  Bessel!  function  of  the  first  kind  of  argument  (k,c/2). 
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FIRST  TERM  IN  THE  ASYMPTOTIC  EXPANSION 
of  V'  AND  V”  As  k,  -  co 


APPENDIX  B 


EVALUATION  OF  THE  ASYMPTOTIC  INTEGRALS 

Section  III  outlined  a  method  for  accelerating  the  numerical  computation  of  Sommerfeld-type 
integrals.  It  was  shown  that  the  computational  efficiency  could  be  considerably  enhanced  by  sub¬ 
tracting  the  slowest  decaying  asymptotic  term  of  the  integrand  and  evaluating  it  in  closed  form. 
The  derivations  of  the  closed  form  solutions  for  the  asymptotic  integrals  are  presented  in  this 
Appendix. 

Consider  the  following  asymptotic  integral,  associated  with  the  matrix  element  Zg?^  (Table  I) 

ou 

A&*  =  I  VA  U*M  Mk.  (Bl) 

m  o 

where 


Um.v  Jm(ktcp)P;(p)dp  (B2) 

*1  I) 

P „'(p)  =  —  jvT^p)  +  (i/+2)T v+2(p)  J/V 1  -p2  (B3) 


-jk« 

W€0(l+€r2) 


(B4) 


"J"  is  a  Bessel  funciton  of  the  first  kind,  and  c  is  the  radius  of  the  circular  patch. 

Substituting  (B2-B4)  into  (Bl)  and  interchaning  the  order  of  integrations  yields  the  following 
expression 


i  i 

A«L  =  C  /  drP/(r)  / dpPj(p) 


/  Jm(k,cp)Jm(k,cf)dkt 

O 


(B5) 


where  C  = 
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O>€0(l+€r2i 
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The  quantity  in  the  !•  •  •}  brackets  is  recognized  as  being  related  to  the  Associated  Legendre 
function  of  the  second  kind  of  integer  half  orders  [6].  or  more  precisely 
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Jm(  ktcp)Jm(k(CT  )dkt  — 


— 7 — \in  ^  1 

ndTp)lu 


p-  +  r  * 


2  pr 


( B6 ) 


The  function  Qm_ l/2  can  be  numerically  evaluated  from  its  recursion  relations  [6]  using  the  well- 
known  back-recursion  algorithm  [14]. 

To  enhance  the  clarity  of  the  presentation  it  is  expedient  to  introduce  the  following  funda¬ 
mental  integral 


i  l 

I(i'.ui.m)  —  J  d T  J dp 

o  o 


Tu(r)TJLp) 

■J  1— r2  V 1  — p2 


p2  +  T2 
2pr 


1  r 

=  f  dTf  dP 


[Tdr)TJp)  +  T„(p)TJt)] 
V 1— r2  J  1—p2 


75? 


Q 


p2+r2 


2  pr 


(B7) 


Utilizing  (B7).  the  expression  for  A^,  can  now  be  rewritten  as  follows 


A  pp  s  C. 
nm.uiv  v'l 


v<j)l(v,(D,m )  +  v(tu+2)I(i',(u+2.m) 


+  (i/+2)wT(i'+2.<d.m)  +  (v+2)(ti>+2)I(v+2.w+2.m) 


where  Cf  = 


~j£ 

7rw€0(l+€r2) 


(B8) 


Related  asymptotic  integrals  have  been  evaluated  in  an  analogous  manner.  The  resulting  expres¬ 
sions  are  summarized  in  Table  III. 

Several  observations  concerning  the  evaluation  of  the  asymptotic  integrals,  and  more 
specifically  the  integral  Ki'.to.m).  should  be  made.  The  most  important  of  these  is  the  fact  that 
I(i'.ciMn)  does  not  depend  on  any  geometrical  or  electrical  parameters  of  the  problem.  Hence  it 
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should  be  computed  only  once  for  each  relevant  Lombination  of  integers  (v.ut.m)  and  then  acred 
for  subsequent  recall.  Furthermore,  utilizing  the  hypergeometric  function  series  representatiea  for 
Qm— 1/2  and  integrating  term  by  term,  the  double  integration  in  (B7)  can  be  reduced  to  an  mCnite 
series  summation. 

The  integral  of  equation  (14)  can  be  evaluated  in  terms  of  elementary  functions.  With 
the  aid  of  Table  II  and  equation  (12).  A™  is  rewritten  as  follows 


A  XX 


-C .// 


oo  jr/2  Sin 

cos^a - 


Aktcosa 


I)  <> 


Aktcosa 


sin^d^wsina) 

(k,wsina)2 


•  cos(Akt|p— q|cosa)  kt2  dktda  =  CxA(|p— q|) 


(B9) 


where  Cx  = 


_A 

IT 


*> 


i 

W€0(€rl+€r2) 


Formulas  on  pages  449-453  in  reference  [6]  can  be  applied  to  evaluate  A(|p-q|).  Considerable  alge¬ 
bra  is  required  to  obtain  the  final  result,  which  is  summarized  in  Table  IV. 

.Note,  evaluation  of  the  asymptotic  integrals  associated  with  Z^^,,  Z^|p  is  omitted  here,  since 
it  was  shown  in  Section  III  that  for  the  case  of  the  electromagneticaliy  coupled  patch  the  integrands 
of  these  matrix  elements  decay  exponentially. 
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TABLE 


m 


Evaluation  of  Integrals  A^,.  A,£^. 


A  PP  =  C. 


vwlivM.m)  +  v(<i>+2)l(x'.w+2.m)  +  ( v+2)u>Kv+2.<D.m )  +  (»'+2)(a>+2)l(i'+2.a»+2.m) 


A,^  =  C2  Ul(w.£.m)  +  (co+2)I(w+2.£.m) 


A***  =  C3l<7).£.m)  +  C< 


I(t)+1.£+1  .m+1 )  +  I(t)+1.£— 1  .m+1 )  +  ICtf—  1  .£+1  .m+1 )  +  K  rj— l.£— 1.  m+1 ) 


Constants: 

c,  =  i-  V'A0  :  C2  = 


t r 


— 2mc  .  .  n  - 
— — —  »  ao  .  v  3  “ 

7 r 


4m  c  v>  ■  n  —  c  \t 

— — —  V  A„  .  I  4  -  —  V  A0 
7T  7 T 


...  _  -j  ....  _ 

V  A,)"  uxJT^)  ■  A°  2 


Notes: 

1)  Kv.ai.m)  =  I(to.v.m)  ;  Kv.— <u.m)  =  Kv.<u,m) 

2)  The  second  term  in  A*^  results  from  subtraction  of  V"A  from  the  original  integrand 
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TABLE  IV 


Evaluation  of  A^*  =  CxA(|p— q|) 


A(0)  -  ^(2K3-K4)  |  +  I  f 

r>  o  24  a=/j  O  os/j 

/j  ^ 

A(l)»  JL(2K4-3K3)  |  +  JL(7k,-21K2+15K3-3K4)  |  + 

24  ar=/0  48  o=/j 

+  -3L(-25k,+27K2-9K3+K4)  |  +  Z-K,  | 

4*4  Or — /j 

7T  ^n~l 

A(n)|n>2=  JL  [(n-2)3K,-3(n-2)2L2+3(n-2)K3-K4]  |  - 

IT  ^ 

-  JL.  [(3n3— 6n2+4)K]— (9n2— 12n)K->+(9n— 6)K3— 3K4] 

48 

TT 

-  [(— 3n3— 6n2+4)kt+(9n2+12n)k2— (9n+6)K3+3k4]  { 


[( n+2  )3K  ,-3(  n+2  )2k2+3( n+2 )k3-k4] 

48  Q 


Definitions 


ii  -  A  1  ..  1  ,  .  a 

1-  I'm  =  — - -  -  ;  k->  =  —  In  tan  —  : 

2w2  sina  ‘  w  2 


»-  _  2  ,  .  tr  ot  4w  1 

K3  =  -  /n  tan  _+_  :  K4  =  — - 

a  4  2  A2  cosa 


2.  ln  =  arctan  :  l0  =  0  ;  =  y 
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GROUND  PLANE 


2  =  0  6  -  --  ■  ■  6 

SHORT  CIRCUIT 


Two-layer  grounded  dielectric  slab  and  its  equivalent  circuit  representation. 


*  FEED  LINE 

\  i  €0'Mo  PATCH 
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Figure  4(b) 
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